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We develop an effective tight-binding Hamiltonian for spin-orbit (SO) interaction 
in bent carbon nanotubes (CNT) for the electrons forming the vr bonds between the 
nearest neighbor atoms. We account for the bend of the CNT and the intrinsic spin- 
orbit interaction which introduce mixing of vr and a bonds between the Pz orbitals 
along the CNT. The effect contributes to the main origin of the SO coupling - the 
folding of the graphene plane into the nanotube. We discuss the bend-related con¬ 
tribution of the SO coupling for resonant single-electron spin and charge transitions 
in a double quantum dot. We report that although the effect of the bend-related 
SO coupling is weak for the energy spectra, it produces a pronounced increase of the 
spin transition rates driven by an external electric field. We find that spin-flipping 
transitions driven by alternate electric fields have usually larger rates when accom¬ 
panied by charge shift from one dot to the other. Spin-flipping transition rates are 
non-monotonic functions of the driving amplitude since they are masked by stronger 
spin-conserving charge transitions. We demonstrate that the fractional resonances ~ 
counterparts of multiphoton transitions for atoms in strong laser fields - occurring in 
electrically controlled nanodevices already at moderate ac amplitudes - can be used 
to maintain the spin-flip transitions. 
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I. INTRODUCTION 

The spins of electrons confined in qnantum dots are considered candidates for the qnan- 
tuni information processing, and the external control of confined spins has attracted a lot of 
attentiorP. In III-V devices the possible applicability is limited by low coherence time due 
to the interaction of the electrons with the nuclear spin fielcP®. The absence of the nuclear 
spin-held in carbon-based materials: graphen^ and carbon nanotubed^ makes them attractive 
for coherent spin manipulation. In contrast to graphen^ in which an electrostatic conhnement 
of charge carriers is excluded by the Klein tunneling, formation of quantum dots by external 
potentials is straightforward in semiconducting carbon nanotubes (CNTs)El. TPe electron spin 
can be controlled with electric helds by the spin-orbit coupling. The spin-orbit coupling in 
graphene is of an atomic origin and couples the pz orbitals (forming the vr bonds) and the 
in-plane Px/Py which form the a between the ions. However, due to the orthogonality of pz 
and in-plane orbitals the spin-orbit interaction does not mix the electron spins in the band 
structur^. This is no longer the case once the graphene plane is folded into the CNT^lSl, for 
which the effect of spin-orbit coupling are clearly observed in spite of low atomic number of 
carborP^. 

The spin control in CNTs is intensely studied by both experimenlP^ and theory®®IISH22] 
Recently, experimental observations of spin transitions driven by external ac electric held were 
reportecpl^ with the electric-dipole spin resonance (EDSR) mechanism that was applied pre¬ 
viously for III-V quantum dotpl^. The effects of the EDSR-driven transition j^^*^^*^^ * ^^*^^ * ^^! are 
observed by lifting the blockade of the current how through a biased double quantum dot. 

The spin-transitions in carbon nanotubes can be triggered by a bend of the CNpH, which 
induces a dependence of the Zeeman splitting on electron position within the nanotube via the 
strong anisotropy of the ehective Lande factor g. A subsequent experimentP^ indicated that the 
bend of the CNT has indeed a dominant contribution to the EDSR ehect. Recent paperP®^] 
investigated the CNT bend in the context of the symmetry of single-electron wave functions in 
external electric heldP*, atomic disordeP', and the leakage current!^. 

In this paper we consider the contribution of the bend of the CNT to the principal spin-orbit 
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interaction resulting from folding the graphene plane into the tub^. We develop an effective 
tight-binding Hamiltonian for the SO coupling resulting from the bend of the CNT as a whole. 
We consider the effects resulting from the bend of the CNT for electron spin manipulation in a 
double quantum dot. The contribution of bend-related SO interaction for the results should be 
expected small since the radii of the bend are by at least an order of magnitude larger than the 
CNT radii. We indeed hnd that the effects of the spin-orbit coupling resulting from the bend 
for the energy spectra are negligible, but the spin transition driven by EDSR are accelerated 
several times. We hnd that generally the spin-hip transitions occur with a larger rate when 
the interdot charge transfer accompanies the inversion of the spin as compared to intradot spin 
hips with hxed charge distribution. The resonant frequencies for transitions involving charge 
transfer with conserved or hipped spin are close also at high magnetic held. We hnd that the 
direct Rabi oscillation with spin inversion becomes attenuated by the spin-conserving interdot 
charge transition for a larger amplitude of the AC held. However, the spin-hip transitions can 
still be reached by the fractional resonances. 

The experimental conditions with ac bias applied to the CNT resemble the ones present for 
atoms and molecules in strong laser helds. In these systems besides the common Rabi oscil¬ 
lation, fractional resonances - also known as higher harmonics generation - are observecP^^, 
which are counterparts of multiphoton transitions of quantum optic^. We indicate that in 
conditions when the transition with spin hip and interdot charge transfer is masked by a strong 
spin-conserving transition, one can employ fractional resonances to perform the spin-hips. We 
indicate that the rates of spin-hipping transitions are non-monotonic functions of the driving 
frequency due to overlap with the strong spin-conserving resonances. The SO coupling due to 
the bend largely increases the width of the fractional resonances for spin-hipping transitions. 


II. THEORY 


We consider a CNT of length L = 31.8 nm - see Fig. [^a). The CNTs geometry is dehned 
by a chiral vectoi® = riiai -|- n 2 a 2 with primitive lattice vectors, ai = ao(l,0) and a 2 = 
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FIG. 1: (a) Schematics of the considered system of a nanotube of length L = 31.8 nm (150 elementary cells), 

radius r = 0.78 nm, with a circular bend of an arc radius i? = 30 nm. Magnetic field is applied in the z direction. 
Along the z direction a double quantum dot is defined and a bias applied, (b) Cross section of the pink ring 
of (a) with the definition of angles 6 used for hopping parameters, and a local coordinate system with Zi along 
the Pz orbital of ith ion. (c) An enlarged fragment of (a) and the coordinate system used to rotate 

the spins defined along the Zj direction, (d) Energy spectrum for B = 0 and bias field Fbias = 10 kV/cm in 
function of the depth of the Gaussian quantum dots. The red and green lines correspond to energy levels which 
are localized mostly in the left and right dots respectively. Each of the energy levels plotted with green and 
red lines is nearly fourfold degenerate with respect to the spin and valley, (e) Energy levels as functions of the 
magnetic field at F^ias = 10 kV/cm. The purple, green and blue arrows show the allowed, i.e. valley conserving 
transitions from the K' f {1) ground state to K' 4 (/), K' f (r) and K' 4 (r) energy levels, where I and r denote 
the localization of the most of the charge in the left or the right quantum dot. (f) Energy levels for V = 0.55 
eV as functions of the bias for B = 0, and charge density for the branches localized in the left and right dot at 
Fbias = 2 kV/cm (g) and Fbias = 10 kV/cm (h). 
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ao(l/2, ■\/3/2) and oq = 0.246 nm. Here we consider a zigzag nanotnbe with rii = 20, n 2 = 0 
for which the CNT radius is r = 0.78 nmP^, and for which an energy gap is present in the 
dispersion relation that allows for confinement of charge carriers by external potentials. We 
assume that the radius of the bend is i? = 30 nm [see Fig. [^a)], unless stated otherwise. 

The atomic spin-orbit interaction is due to the electric field of the nuclei. The spin-orbit 
interaction due to this field mixes the pz orbitals with p^, and Py orbitals of opposite spins. The 
mixing leads to hybridization of pz orbitals [labeled as zj for j-th ion - see Fig. 1(a)]- forming the 
71 bonds with the p^/y (labeled as Xj,yj) - forming the a bonds between the nearest neighbor 
atom^^ ^^Fsl ^ For | I izj) standing for the spin-orbitals along the Zj orbital of jth ion 
the first-order corrections to the spin-orbital \zj ^zj) due to the atomic spin-orbit interaction 
stemming from the electric field of the nuclei have the fornP^*!^ 

\^j Tzj) ~ \^j Tzj) S\Xj -\- ipj -l-zj) (1) 

and 

\^j '['Zj) ~ \^j '['Zj) T S\Xj ipj Tzj); (2) 

where 6 parameter describes the strength of the SO coupling. For the fiat graphene the and 
x/y orbitals for the nearest neighbor atoms are orthogonal, but the orthogonality is lifted when 
the graphene plane is folded into the CNT^SHso], introduces spin-orbit coupling effects 

to the 71 bancP. 

The second-order tight-binding Hamiltonian theories taking into account the 2s,2px,2py, 
and 2pz orbitals are availabl d^^*^^* ^. The contribution of the entire basis introduces combined 
effects including the electron-hole asymmetry in the spectrumP^® and spin-splitting of the 
bands for chiral CNT’d^^*^. In this paper we restrict the study to zigzag (achiral) CNT and 
consider the conduction band states only. We aim at estimation of the spin-orbit coupling 
effects due to the bend of the entire CNT. Since the effects of the bend should be expected 
much smaller than the ones due to folding of the graphene to the CNT we restrict the modeling 
to the lowest-order corrections only®2]_ Namely, we consider only the pz electrons, and the 
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effects of mixing of type given by Eqs. 0 and (|^. 

We look for an effectiv#*^ tight-binding Hamiltonian of form 

»= E (4, • iy'• + A.C,) 

+ + (3) 

i,cr,cr' 

The hrst summation runs over spin-orbitals of nearest neighbor atoms, (qo-) creates 
(annihilates) the electron at ion i with spin a in z direction, and is the spin-dependent 
hopping parameter. The second summation in Eq. ([^ accounts for the external potential and 
the Zeeman interaction with the Lande factor g = 2. cr stands for the vector of Pauli matrices, 
and B for the magnetic held vector. The magnetic held is necessary for the EDSR transitions 
to be observed. We assume that B is applied in the y direction [see Fig. 1(a)]. 

The double quantum dot (DQD) potential induced by e.g. external gates is modeled as a 
sum of two Gaussians Vj:,qd{z) = —V [exp((y — sY/dP) + exp((y -|- sY/d?)], where the distance 
between the centers of the dots is 2s = 10 nm, and the width of a single QD is 2d = 4.4 nm. 
Note, that the DQD potential - depends on the global y direction [see Fig. [^a)] and not on a 
coordinate along the length of the bent CNT. Thus, the bend lowers the angular symmetry of 
the eigenstates with respect to the axis of the tube, which as discussed in Ref.^ allows for the 
spin-flip transitions in presence of the SO coupling. The external potential dehned within the 
CNT in the absence of the ac held is then W = Vogoiz) + eFhiasZ, where the second term is 
responsible for the bias held applied to the system. 

The previous workP^ used the hopping parameters for a straight CNT, while the present 
paper deals with ehects of the bend to the parameters including the SO interaction. To include 
the ehect of the bend into SO coupling we need to modify the hopping parameters. For 
necessary transformations we use four diherent coordinate systems (see Fig. 0 a-c)): (xyz) - 
global coordinate system, {xji/jZj) - local coordinate system with zj aligned with spin-orbital 
at j-th ion and yj parallel to local nanotube axis, (x^y^zY) - auxiliary local coordinate system 
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with coinciding with the global y axis and coinciding with yj axis, [x'^y'^z'^) - system 
equivalent to XjyjZj but rotated by 0 angle about the x^ axis. 

The procedure for derivation of the hopping parameters for the bent nanotube is the follow¬ 
ing; 

i) We consider each of the cross section of the nanotube (pink line on the Fig. [^a)), separately 
and start from the eigenspinors | tzj), I Xzj) defined along the zj orbital - as in Ref.® [see Fig. 


ib)|. 


ii) We rotate these eigenspinors by the 0 angle about the axis = (cos6*j, 0, sin^j) (see Fig. 
Sb- c)). The rotation matrix has the form; 




cos ^ — i sin ^ sin 6j 


—i sin ^ cos 6j 


-i sin % cos 6j cos ir + i sin % sin 6j 


(4) 


with (Tj - the Pauli matrices in the local coordinate system xjyjZj. We get new eigenspinors 
in the rotated system x'^y'^z'^ with axis y' coinciding with the global direction; 


I ) = (cos ^ - z sin ^ sin Oj)\ ) - z sin ^ cos Oj\ ) (5) 

I ) = -z sin ^ cos 9j \ (cos ^ z sin ^ sin Oj) \ 4 .^.) (6) 

in) We transform the eigenspinors | I 'i'^O along the z', direction in the xW,z', coordinate 
system into eigenspinors in z direction in the xyz coordinate systemp! | '(z), | [see Fig. 
[3a,c)]; 

pi6jl2 

I t^) = + *1 izp} 

= ^{(cos|e*^^/^ + sin|e-*^^/^)|t.,) 

+ z(cos 


(7) 





iv) Using transformations in) we convert formulae Q and (|^ to the \zj tz) basis: 


lU't^) ~ \zj t^) + iS cos(()j\xj tz) - iS sm(()j\xj 1^) 

+ i6 sin (f)j sin 9j\yj (9) 

+ ^(cos^ + sin^ 

and 


\zj iz) ~ \zj iz) — sin (f)j\xj tz) — cos (j)j\xj J,^) 

- (5(cos2 + sin^ t^) (10) 

— i6 sin (f)j sin 9j\yj 

v) Finally, using iv) and discarding higher-order terms in 6 we obtain the hopping parameters 
to be used in Hamiltonian (1); 

= {-.i t. \H\^i t.) 

= {zi\H\zj) + i5 cos (l)j{zi\H\xj) — i5 cos (l)i{xi\H\zj) 

+ i5 sin (l)j sin 9j{zi\H\yj) 

- i5sin(i)iSm.9i{yi\H\zj), 


( 11 ) 
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ti} = {zi U \H\zj U) 

= {zi\H\zj) — i6 cos (l)j{zi\H\xj) + i6 cos (l)i{xi\H\zj) 

— i6 sin (pj sin 9j{zi\H\yj) 

+ i6 sin (pi sin 9i{yi\H\zj), (12) 


and 


= {Zi t. \H\^j U) 

= — i6 sin (j)j{zi\H\xj) + i6sin (()i{xi\H\zj) 
- 6{sin^ + cos^ Ye-“'‘){zi\H\yj) 

+ i(sin^ ye* + cos^ ye“*)(!/il-ff|2j), 


(13) 


(j; = {Zi U \H\Zj t4 

= —i5 sin (j)j{zi\H\xj) + i5sin (pi[xi\H\zj) 

+ 5(sin2 + cos^ 

- + cos^ ^e^^^){yi\H\zj). (14) 

Following RefP the matrix elements for the neighbor pz orbitals to be used in the hopping 
parameters above read 


{ai\H\aj) 


= K>(«*) ■ n(aj) 


+ [V^ - V^) 

' ^ pp 


(n(ai) ■ %)(n(aj) ■ %) 


ppj 


IR 


ji\ 


(15) 


where a = x, y oi z, ai is orbital localized at site R,; (Fth ion), 11 ( 0 ;,) is a unit vector in the 
direction of orbital a* [see Fig. I3b)l- 

For a straight CNT, (pj = 0 for all j, the formulae reduce to the form given by Ref.^, for 
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which non-zero spin-flipping hopping parameters appear only along the circnmference of the 
nanotube {6i 7 ^ 6j), where the Pz orbitals are no longer strictly parallel. The form of spin-orbit 
coupling derived above accounts for the deflection of Pz orbitals along the length of the bent 
CNT, and the spin-flipping hopping parameters appear also for neighbor orbitals with same 6 
but varied 0. In the numerical results below we use the spin-orbit coupling parameter 6 = 0.003 
after Refs .121121. and the tight-binding parameters = —2.66 eV and = 6.38 of Ref.^ 

The magnetic field induces the spin Zeeman effect in Hamiltonian (1). The orbital ef¬ 
fects of the magnetic field are introduced by Peierls phase to the hopping terms tij{B) = 
f,j(0)exp(ig//^^ A-dl), where B = V x A, and = h/e is the flux quantum. 

Below we discuss the contribution of the SO coupling due to the bend. The reference results 
are obtained for the hopping parameters as for the straight CNT® and the bend of the tube 
enters only the potential energy and the vector potential through the spatial coordinates of the 
ions. 


III. STATIONARY STATES 

The calculated energy spectrum of the CNT is plotted in Figure [^d) for Fbias = 10 kV/cm. 
With the black points we plotted the energy levels of the states which are localized outside 
the double quantum dot potential as functions of the depth of the Gaussian cavities. The two 
horizontal black lines correspond to energy levels localized at the zigzag edges of the tube. 
The color lines indicate the energy levels localized within the DQD, the red (green) ones with 
the majority of the charge localized in the left z < 0 (right .^ > 0) dot. Below we assume 
V = 0.55 eV, for which one red and one green energy levels - each nearly four-fold degenerate 
with respect to the valley and the spin, appear below the neutrality point energy (zero). We 
will consider a single-electron localized in DQD inside a CNT which is otherwise neutral. 

In Figure [^f) we plotted the energy levels as a function of the bias field. Each of the plotted 
energy levels is two-fold degenerate. At F^as = 0 the splitting between the green and red 
energy levels is due to the tunnel coupling between the dots and formation of bonding and 
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antibonding-orbitals. Each pair of the energy levels is split by the spin-orbit coupling energy. 
As the bias is applied the lower (upper) quadruple of energy levels becomes localized in the left 
(right) dot. The electron localization of the energy levels is given in Figure [^g-h) for the bias 
helds Fbias = 2 kV/cm and 10 kV/cm. 

The magnetic held splits the degeneracy of the energy levels - see the energy spectrum of 
Fig. 13 e). We label the states by the valle}^the majority spin component and the dot left 
or right {l,r) in which most of the charge is localized. Below we consider transitions from the 
ground-state K' j" to the excited states at i? = 5T. The present calculation accounts for both 
the valleys, but since the modelled CNT is defect-free no intervalley transitions are observed. 
We discuss transitions from K' ground-state to the hnal states of the K' valley. 

The spin-orbit coupling which stems from the hnite radius r of the nanotube results in 
formation of the spin-valley doublets split by the spin-orbit interaction energy A^o- The 
energy effect of the spin-orbit coupling resulting from a hnite radius of the bend R is too small 
to be plotted in Fig. [^e) and for i? = 5T it is of the order of 0.01 meV at most, depending 
on the energy level. A relative effect on the average spin of the states is larger. For instance 
{sz) for the ground state at 5T is reduced by the spin-orbit coupling due to the bend from 
0.4991h to 0.4805h. Figure compares the ground-state spin and charge density obtained for 
SO coupling for a straight CNT with the ones including the effect of the bend. The charge and 
spin distribution are similar in shape, but the minority spin density is increased by about 25 
times by the SO due to the bend. 


IV. SPIN AND CHARGE TRANSITIONS 


We investigate the spin and charge dynamics of the system when an external ac electric 
held is applied along the 2 ; direction, VAcif) = eFAcz?Fa-{ojt). For that purpose we solve the 
Schrodinger equation ih^ = hf'T for Hamiltonian H'(t) = H + VAc{t) in the basis of H 
eigenstates, 

N 

T(r,(T,f) = ^Cn(t)T„(r,(T)e 


n=l 


( 16 ) 
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FIG. 2: Charge (c,f) and spin (a-e) densities for the two-electron ground state at Fi,ias = 10 kV/cm and 
B = 5T. Results correspond to the bend nanotube. In (a-c) SO coupling (hopping parameters) was adopted 
for a straight CNT. The SO that accounts for the bend was used in (d-f). Arbitrary units for the probability 
density are used. 


where H^n = En^n- For the wave function (16) introduced to the Schrodinger equation, after 
application of the Galerkin projection, the time evolution is given by a system of ordinary 


differential equations 


N 

Cn{t)eFAc sm{ut){^k\z\'^n)e « , (17) 

n=l 


ioT k = 1,2,... N that we solve with the Crank-Nicolson time-stepping method. We consider 
that all the energy levels below the ground-state level localized within the DQD, with energies 
E < —300 meV for V = 0.55 eV [see Fig. 1(d)] are fully occupied. For the time evolution we 


take N = 12 lowest-energy spin-orbitals to the basis (16) starting from the DQD ground-state. 


The results can be considered as exact solutions of the time dependent Schrodinger equation 
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(a) ^K'T(r)^K'i(r) (c) ^K'l(/) ^ K'T (r) ^ K'i (r) 
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(e) ^K'HD ^K’T(r)-K’i(r) 
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(g) ^K'i(Z) ^K'T(r)^K'i(r) 





FIG. 3: Upper panels of (a-h): maximal occupation probability for the three K' excited energy levels reached 
during 500 ps time evolution from the K' f (/) ground state in the initial condition. Pink/green/blue lines 
correspond to the transitions to K' {l)/K' f {r)lK' (r) states, respectively. The triangles on top of the 
upper axis indicate the energy difference between the ground-state and the excited K' states as calculated from 
the energy spectrum. The lower panels of (a-h): minimal spin obtained during the evolution. Light green lines 
indicate minimal spin for 500 ps simulations (corresponding to results from upper panels), red/black lines show 
intermediate results for 200 ps/100 ps simulations. The applied bias is Fbias = 2 kV/cm (a-b,e-f) and Fbiae = 10 
kV/cm (c-d,g-h). The amplitude of the ac field is Fac = 0.5 kV/cm in (a-d) and Fac = 4 kV/cm in (e-h). Plots 
(a,c,e,g) correspond to the results obtained with the hopping parameters neglecting the bend and (b,d,f,h) with 
hopping parameters including the bend. The plotted scans were calculated with a hw spacing of O.lfieY. 
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since introduction of a larger number of eigenstates to the basis does not alter the results in a 
detectable manner. 


A. Spin-conserving transitions 

The ac electric field drives the transitions between the four K' valley states. Figure]^ shows 
[upper panels of Fig. |^a-h)] the maximal projection (|(4/(f)|\['„)p) of the evolving quantum 
state starting from the ground K' j' (Z) state on the three excited K' states for simulation that 
covers 500 ps. The green line in the upper panels of Fig. [^corresponds to the spin-conserving 
transition K' (1) —>■ iF' j' (r) with the electron passing from the left to the right quantum dot. 
The spin-conserving charge transition is very fast. At the bias voltage Fbias = 2 kV/cm for the 
resonant energy hu equal to the energy difference between the K' j' (/) and K' j' (r) energy 
levels the charge transition from the left to the right dot appears within about 8 ps {F^c = 0.5 
kV/cm). For Fac increased to 4 kV/cm the corresponding transition time is about 1 ps only. 
The charge transition between the dots slows down when the bias increases the detuning of 
the single-dot energy levels and reduces the overlap of the spatial wave functions of the initial 
and hnal states. For F^ias = 10 kV/cm the charge transition times are ~ 14 ps and ~ 2 ps for 
Fac = 0.5 kV/cm and Fac = 4 kV/cm, respectively. The lower rate of the charge transition 
is consistent with the dipole matrix elements of Table I, and is found also for the CNT with 
SO including the effect of the bend. Note, that the amplitude of the ac held has a pronounced 
inhuence on the width of transition K' 'I (1) ^ K' 'I (r) [cf. Fig. [^a,e) and Fig. ic ,g)]. For 
a small amplitude of the alternate electric helds Fac = 0.5 kV/cm we observe [upper panel of 
Fig. |3ta.c )] a narrow resonant transition at half the energy spacing between the ground state 
K' j' (Z) and the excited state K' (r) (see the peak marked by 1/2 in the Figures). This is 
the fractional resonance, which is a counterpart of two-photon transitions that are observed for 
atoms and molecules in strong laser held^. When the amplitude of the ac held is increased, one 
observes a series of higher-order resonances with the resonant amplitude up to 6th harmonics 
[upper panel of Fig. |^e, g)l- 


15 


R [nm] 

Fbias 

K' i (/) 

K 't (r) 

K' i (r) 

15 

2 kV/cm 
10 kV/cm 

0.068 (0.347) 
0.016 (0.074) 

4.815 (4.666) 
2.852 (2.340) 

0.180 (1.078) 
0.317 (1.565) 

30 

2 kV/cm 
10 kV/cm 

0.032 (0.164) 
0.009 (0.043) 

4.847 (4.820) 
3.010 (2.877) 

0.090 (0.469) 
0.175 (0.859) 

120 

2 kV/cm 
10 kV/cm 

0.008 (0.040) 
0.002 (0.011) 

4.857 (4.856) 
3.059 (3.051) 

0.023 (0.112) 
0.044 (0.220) 


TABLE I: Dipole matrix elements \ {K' f (/)|ez|T)| for the transition from the ground state to one of the three 
excited states of the same valley, in units of ex nm for S = 5 T. Results without (with) the inclusion of the SO 
coupling due to the bend of the nanotube with radius R are given outside (inside) parentheses. 



E [meV] 


E [meV] 


FIG. 4: Spin conserving (black lines) and spin flipping (red lines) transition probabilities for a varied duration 
of the AC field. In (a) the hopping parameters of a straight CNT are used. The bend of the CNT is included 
in (b). We use Fuas = 10 kV/cm and Fac = 4 kV/cm as in Fig. 3(g,h). 


B. Spin-flipping transitions 


Table I lists the dipole matrix elements for the transitions, which are inversely proportional to 
the transition times. The oscillator strength for the transition with conserved spin K' ^ (r) 
exceeds by orders of magnitude the ones for spin-flipping transitions. For a larger bias an 
increase of the matrix element for transition with spin-flip accompanied by interdot charge 
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hopping K' (1) ^ K' I (r) is observed. An opposite tendency is found for the intradot spin 
flip K' ^ (/) —)■ A'' 4 , (/). In consistence with the data for the matrix elements given in Table 
I, the time-dependent calculations indicate that for low amplitude of F^c = 0.5 kV/cm [Fig. 
[^a,c)] the interdot spin flip is by far a more pronounced transition than the spin flip with the 
electron staying in the left dot. For SO coupling taken for the straight CNT the latter is weak 
at low bias [Fig. |^a)] and nearly absent for the larger bias [Fig. |^c)]. 

The fact that the spin-flip interdot transition occurs faster then the intra-dot spin inversion 
is quite counterintuitive, but can be explained in simple terms based on the approximate sym¬ 
metries of the wave functions. We hnd that the lower-energy states that are mostly localized in 
the left dot have a bonding character, i.e. both the majority t and the minority j, components 
have a positive average value of the z-parity operator dehned as Pzf{x, y, z) = /(x, y, —z). 
On the other hand, the K' j'J, (r) excited states have an antibonding character of both the 
spin components with a negative average value for P^,. The transitions matrix elements (Table 
I) between pairs of states of opposite character (bonding-antibonding) in terms of formation 
of artihcial molecular orbitals are naturally larger than for the states of the same (bonding¬ 
bonding) character. Therefore, the intradot transition K' ^ (/) —)■ iF' j, (r) has a larger matrix 
element than the interdot one K' ^ {1) ^ K' {1). 

Let us now compare the results as obtained for SO interaction as set for a straight CNT 
|Fig. |3|;a.c,e ,g)], and the ones for the SO theory accounting for the bend [Fig. 3(b,d,f,h)]. 
The spin-orbit interaction due to the bend has no significant influence on the rate or width 
of spin-conserving transitions. The bend-related SO interaction shortens the direct transition 
rate to the K' j, state - see the minimal (s^) for 100, 200 and 500 ps simulation. For Fuas = 10 
kV/cm the transition to K' j, (r) state is displayed in a closer detail in Fig. as a function 
of time for SO as obtained for a straight and bent CNT. The transition time is reduced 5 
times by the contribution of the bend to the SO interaction, in consistence with the results for 
the matrix elements given in Table I. In Fig. |^we plotted the induced transition probability 
without (a) and with (b) the bend-related contribution to SO interaction for parameters of Fig. 
3(g,h) as a function of time. For identihcation of the separate lines - see Fig. 3(g,h). The 
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fully developed spin-flip and charge hopping line that is related to the half-resonant transition 
(see the thin red line for E = 6.3 nieV) requires 20 ps [Fig. |^a)], while the theory neglecting 
the bend contribution to the SO interaction predicts transition within 100 ps [Fig. |^b)]. The 
intradot spin-flip transition line (the red peak near 0.4 eV in Fig. |^a,b)) as well as fractional 
(1/3) resonant transition with spin flip and the charge hop (the red peak near 4.2 nieV) is 
fully developed within 100 ps (full theory) time versus 0.5 ns (neglected SO contribution due 
to the bend). The matrix elements for the spin-flipping transitions increase with account taken 
for SO effects due to the bend since the contribution of the minority spin component to both 
the initial and the hnal states are made larger by the bend-related SO contribution [Fig. |^. 
The increased coupling between the states of opposite spin orientation results also in the larger 
width of the spin-flipping transitions [cf. Fig. [^. 

Table I presents the results also for the radii of the bend i? = 15 nm and R = 120 nm. The 
spin-conserving transition rate is independent of the bend {R). The spin-flipping rates decrease 
with R - with or without inclusion of the bend to the SO interaction. For a straight CNT 
the spin-flipping transitions do not occur due to the angular orthogonality of the initial and 
hnal states unless an effect lowering the rotational symmetry is presenlP^ (an atomic defect, a 
perpendicular electric held etc.). Independent of R we hnd a 5-fold increase of the transition 
rate by the contribution of the bend to the SO coupling. 



FIG. 5: K 't (1) K' (r) transition: square of the absolute value of the wave function projection on the 
K' 1 (r) state as a function of time for bias F = 10 kV/cm and the driving amplitude Fac = 0.5 kV/cm. 
The dashed (solid) curve corresponds to the results obtained with the hopping parameters including the bend 
(neglecting the bend). The frequency was tuned to resonance - see Fig. [^c-d) for each case. 
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C. Strongly driven system 

For the larger amplitude of the AC held, F^c = 4 kV/cm [Fig. e-h)], the evolution of the 
system becomes non-perturbative and the matrix elements given in Table I loose their direct 
relation to the actual dynamics of the quantum system. In particular, the direct - hrst order 
- transition of the Rabi type K' 'j' (Z) K' I (r) is missing in the Figure. The width of the 
spin-conserving transition K' (1) K' (r) is largely increased with the amplitude of the 
AC held. This spin-conserving and the spin-hipping transition K' 'j' (Z) —>■ iC' (r) appear close 
to one another on the energy scale. For large Fac the spin-conserving resonance extends over 
the frequency characteristic to the spin-hipping one. Then, a competition between the hnal 
transition states appear and the one which prevails is the spin-conserving transition of a much 
larger oscillator strength [Table I]. For Fuas = 10 kV/cm and Fac = 4 kV/cm the spin-hipping 
transitions for fvjj = 12 meV are only observed for the SO including the bend of the nanotube 
[Fig. [^h)]. Note, that although the direct transition for the nominally resonant frequency to 
the right dot with the spin-hip is absent for larger Fi^i^s = 2 kV/cm and Fac = 4 kV/cm, the 
hnal state K' (r) can be reached with the second- (1/2) and third- (1/3) order transitions 
as given in the Fig. [^e-f). For Fac = 4 kV/cm the half-resonant peaks for the transitions to 
K' J, (r) are outside the half-resonant peak for the K' (r) transition. 

Note, that the direct transition K' '[^ (Z) —?• iF' (r) appears faster for a larger bias - see 
the lower panels to Fig. [^a-d) which show the minimal spin observed for simulation lasting 
100, 200 and 500 ps. This is on the contrary to the behavior found for the spin-conserving 
K' t (Z) —)■ iF't (^") transition which slows down for a larger bias. 

In Figure [^ we plotted the maximal probability to hnd the electron in spin down states 
obtained during 50 ps time evolution in function of the amplitude of the AC held. For each 
transition the frequency was set to resonance. We note, that i) for the spin-hipping transition 
within the left dot (—)■ K' (Z)) the transition rate is a monotonic slowly growing function of 
the AC held amplitude; ii) the spin-hipping transitions accompanied by charge transfer from 
the left to the right dot is generally much faster and in) these transitions are non-monotonic 
function of the amplitude of the AC held - similar behaviour has been shown for DQDs in Ref. 
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FIG. 6: The rate of spin-flipping transitions: the maximal |('I'(t)|i4r' 1 (0)P 1 (?"))P obtained for 

50 ps time evolution for initial state ^*(0) = iti' f (?) as a function of Fac- For each value of Fac the frequency 
w was tuned to the resonance. A bias F^ias = 10 kV/cm was applied. Direct spin-flip transition within the left 
dot is plotted with the green points. The other points correspond to spin-flip transition with charge transfer 
to K' 1 (r) state realized in the direct (black color), and fractional resonances (1/2 - red color, and 1/3 - blue 
color). In (a) the SO interaction for a straight CNT is applied. The SO interaction that accounts for the bend 
is used in (b). 

[38l Here, the reason for the non-monotonic behavior is the neighborhood of the spin-conserving 
transition in the frequency domain. As Fac increases the spin-conserving transition; both the 
direct one and the fractional ones increase radically in width. Once, the spin-flipping transition 
hnds itself within the wide peak of the spin-conserving transition, the spin is no longer flipped 
as the AC held is applied; only the charge oscillation between the dots is observed. The rate 
of the half-resonant transition to K' (r) state decreases for Fac > 5.5 kV/cm, when it is 
consumed by the half-resonant spin-conserving transition. Then a higher transition rate can 
be obtained for the 1/3 resonance etc. For the bend-related SO coupling present the maxima 
of the probabilites reach unity, and for lower values of Fac- 
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D. Transitions for the electric potentials applied along the CNT direction 

The results presented so far were obtained for the DQD potential defined along the ^ direc¬ 
tion. Let us consider - a less realistic case - when the DQD potential and bias field are dehned 
along the axis of the bent CNT, by e.g. bent gates that follow the shape of the CNT. The 
results for the transition matrix elements are given in Table II. The elements for spin-flipping 
transitions with SO as taken for a straight CNT - are drastically reduced as compared to Table 
I. For these values the only effect triggering the spin-flip transition is the variation of the angle 
between the B vector with respect to the local axis of the CNT. The bend-related SO contribu¬ 
tion (the figures in brackets in Table II)- breaking the rotational symmetry of the Hamiltonian 
- increases the spin-flip transition rates by as much as 40 times. The spin-conserving transition 
remains insensitive to the way the external electric field is introduced. 


R [nm] 

Fbias 

K' i (/) 

K 't (r) 

K' i (r) 

30 

2 kV/cm 
10 kV/cm 

0.003 (0.125) 
0.0008 (0.035) 

4.858 (4.841) 
3.062 (2.977) 

0.009 (0.357) 
0.017 (0.682) 


TABLE II: Same as Table I, only for the double dot potential and the bias field defined along the axis of the 
bend CNT and not along the global z direction. In each cell of the column we provide two figures. The one 
outside (inside) the parentheses corresponds to the SO coupling introduced for a straight (bent) CNT. 


E. CNT with a defect 

The results presented above were obtained for a clean CNT with a perfect crystal structure 
for which the inter-valley transitions are absent. In order to estimate the valley mixing effects 
for a CNT with crystal defects we have taken the parameters of hgure |^f) and removed one 
carbon ion at a distance of 2.5 nm from the left edge of the CNT. The results are displayed in 
figure The rates of all valley-conserving transitions within K' remain almost unaffected by 
the presence of the defect. The inter-valley transitions which are forbidden for a clean CNT and 
which activated by the defect produce a sequence of narrow peaks in the upper panel of figure 
Already in hgure [^f) one could notice that the the peaks which corresponded to both charge 
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and spin transitions K' "I' (Z) K' (r) were accompanied by peaks-at the same frequency- 
for the intradot spin flips K' 'I' (/) K' (1). Similar effect is observed for the intra-valley 
transitions-which occur at same frequencies for charge conserving and charge hopping processes 
(see the orange and gray lines in the upper panel of hgure [^. As far as the spin transitions 
are concerned the presence of the defect induces an appearance of three additional narrow lines 
with respect to hgure [^f)-one at hen = 4.7 meV with the half-resonant transition to K I (1) 
(inverted spin, valley and dot, light green line in hgure and another one to K I (/) at left 
of the strong transition to K' (Z) near hen = 1.4 meV and half-resonant one near hen = 0.7 
meV. These two transitions occur at relatively low rates: the lower panels of hguresj^f) and[^ 
remain nearly identical up to 100 ps, only after the three narrow spin-transition lines evolve in 
the spectrum. Concluding, for CNT with a defect the intervalley transitions are observed. The 
transitions are narrow at the driving frequency scale and produce only a slight modiheation to 
the spin hip dependence. 



0 2 4 6 8 10 

htn [meV] 


FIG. 7: Same as figure |^f) but with one carbon atom missing at a distance of 2.5 nm from the left edge of 
the CNT. 
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V. SUMMARY AND CONCLUSION 

We have developed the tight-binding Hamiltonian for a carbon nanotube accounting for the 
spin-orbit coupling due to a-n hybridization of the covalent bonds resulting from both folding 
the graphene sheet to a nanotube and the bend of the nanotube. For discussion of the effects 
of the SO coupling contribution resulting from the bend we considered the electron transitions 
within a double quantum dot defined in a clean carbon nanotube in AC electric fields. We 
demonstrated that although the bend-related SO interaction has a negligible influence on the 
energy spectra its impact on the spin-flipping transition times is pronounced. We found that 
for a large amplitude of the AC voltage the spin-conserving transition peak evolves into a wide 
maximum consuming the interdot Rabi oscillation for the spin-flipping transition. We discussed 
the fractional resonances - the solid state counterparts of the multiphoton transitions. We 
demonstrated that for higher AC amplitude the fractional resonances can be used to perform 
the spin-flips with interdot charge transfer. 
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